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Abstract 

Modelling of the point spread function of the UKIRT IRCAM3 array was conducted 
in order to test for any extended emission around the X-ray binary Cyg X-3. We 
found that the point spread function cannot be represented by a simple Gaussian, 
but modelling of the stars required additional functions, namely Lorentzian and 
exponential components. After modelling for the PSF, we found that Cyg X-3 could 
be represented by two stellar-type profiles, 0.56" apart. 



1 Introduction 

Cyg X-3 has been studied for 30 years after its discovery in 1966. During that 
time, data from the radio, infrared and X-ray wavelengths has shown it to 
be a highly unusual source. This paper concentrates on the deepest infrared 
image of Cyg X-3 (Fender & Bell Burnell, 1996) and subsequent modelling of 
the profile of Cyg X-3. 

The images were taken by Fender using the IRCAM3 array at the UKIRT 
telescope in July 1994. IRCAM3 is a 256 x 256 InSb imaging array, operating 
in the 1-5 /im wavelength range. The data here are composite images from 
the four bands: H, K, L and J. 
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2 H,K,L, J-band images 



Figure 1 shows the four bands, H, K, L and J in the wavelengths 2.2, 1.65, 
1.25 and 3.45 /mi respectively. The images represent 355 individual frames 
with integration times between 2 and 25 seconds, binned up for a total inte- 
gration time of 800 seconds. Stars in the K-band have fluxes ranging from 12 
magnitude (Cyg X-3) to 18th mag for some of the faintest stars. The positions 
of the stars used in this paper were from north to south: Cyg X-3, Star Z and 
Star D. These are the three brightest stars in the images. 
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2.1 k-band image 



The k-band image, figure 2 or the finder chart provided by Fender (Fender & 
Bell Burnell 1996) was used for all the modelling done for the IRCAM PSF 
and the IR emission we are presenting. The points of interest are a slight 
extension to the west of Cyg X-e and that the image is not quite circular. 



k— band 




pixels 



Fig. 2. K-band image 

Modelling was then done using the Starlink PISA programme. A typical stellar 
profile from this image was obtained using the southern star, D. It is impor- 
tant to use a star which is well separated from the rest of the field as we are 
important in the faint extended emission. This template star has to also be 
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Parameter Value 



r 0.32 

a 1.89 

Q 0.066 

R c 2.02 
Table 1 

Parameters fitting the profile shown in figure 3. 

bright enough so its image covered enough pixels. Star D was a good compro- 
mise between brightness and separation. With that as a template PISA fitted 
a profile as in figure 3. Parameters to the profile are shown in table 1, the 
meaning of which are explained below. 
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Fig. 3. Profile fitted to trS^rfc^Pfr^^ 'D' in the K-band image. 

These parameters correspond to a set of equations shown below. You see that 
the profile fits very well to the data. RA wobble in the telescope is shown 
in the profile plot by every 3rd point being slightly brighter that the mean 
profile. 



3 Profile equations 

The equation governing the profile and also the PSF of the IRCAM3 array 
can be expressed as: 
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For r > R c 



where 

r - Percentage of peak for crossover where the function changes from Gaussian 

to exponential. 
a - Gaussian sigma. 

Q - Percentage of Lorentzian component. 

i 

R c = a (In (1/r)) 2 - Changeover radius from Gaussian to exponential. 
The PSF can be modelled by three parts: 

> Central Gaussian + Lorentzian modification 

> Exponential tail-off + Lorentzian modification 

inside a given radius, R c , (2 pixels) the profile is modelled by a Gaussian with 
the Gaussian being modified by Lorentzian wings. Outside a radius of 2 pixels, 
the profile has an exponential tail-off with the same Lorentzian wings. 

Using the profile, PISA then searched the k-band image for all objects that 
had profiles of the template type and subtracted them from the image. After 
subtracting, PISA then re-fitted the image for any underlying stellar images 
and figure 4 shows the results of this. 

Here PISA found emission out to a lower level (represented by the size of the 
circles) for all of the stars identified on the finder chart, except the star to the 
far west. We think that is because of the small dynamic-range for that star 
and even thought the peak is similar to other stars found, the background is 
higher and so it was not found. 

The most interesting thing about this fit is that after re-iterating the image, 
PISA could fit an additional component to the Cyg X3 image, either in fore- 
ground or background. We are not sure whether this is a star co-incidental 
with Cyg X3 or whether it is an artifact of PISA. 
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Fig. 4. pisa fit for all 'stellar' images to the k-band image using the PSF described 
above. Units are pixels, 1 pixel=0.286" 

4 Results and conclusions 



This paper attempts to model the UKIRT IRCAM3 point spread function 
based on the finder chart and H,L,J-band images. We find that a simple 
Gaussian function is not sufficient in modelling the PSF, and that one needs 
to include other functions. We find that a composite function of Gaussian 
with exponential tail-off together with a Lorentzian modification is sufficient 
in modelling the PSF. 

When the PSF is applied to the K-band finder chart, one can fit two separate 
images to the image of Cyg X-3. These two images are separated by 0.56" and 
have a flux ratio of 11:1 (left — right). 
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